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ABSTRACT 

Using the first three and a half years of observations from the Energetic Gamma 
Ray Telescope (EGRET) on board the Compton Gamma Ray Observatory (CGRO), 
phase-resolved analyses are performed on the emission from the three brightest high- 
energy 7-ray pulsars, Crab, Geminga, and Vela. For each pulsar, it is found that there is 
detectable high-energy 7-ray emission above the galactic diffuse background throughout 
much of the pulsar rotation cycle. A hardness ratio is introduced to characterize the 
evolution of the spectral index as a function of pulsar phase. While the hardest emission 
from the Crab and Vela pulsars comes from the bridge region between the two 7-ray 
peaks, the hardest emission from Geminga corresponds to the second 7-ray peak. For 
all three pulsars, phase-resolved spectra of the pulse profile components reveal that 
although there is a large variation in the spectral index over the pulsar phase interval, 
the high-energy spectral turnover, if any, occurs at roughly the same energy in each 
component. The high-energy 7-ray emission from the Crab complex appears to include 
an unpulsed ultra-soft component of spectral index ~ —4.3 which dominates the total 
emission below 100 MeV. This component is consistent with the expected emission 
from the tail end of the Crab nebula synchrotron emission. 

Subject headings: gamma rays: observations — pulsars: individual (Crab, Geminga, 
Vela) 
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1. Introduction 

High-energy 7-ray pulsars have proven to have 
complex, often extended, pulse profiles. By ex- 
tracting only those photons which arrive dur- 
ing a specified phase interval, the emission char- 
acteristics of an individual pulse pulse profile 
component can be determined, providing there 
are sufficient photon statistics. Fortunately, the 
thiiee brightest cumpacl sources in the 7=ray sky 
are the high-energy 7-ray pulsars Crab (PSR 
B0531+21), Geminga (PSR J0633+1746), and 
Vela (PSR B0833-45). The spectral proper- 
ties of the pulse profile components from these 
pulsars have been investigated previously us- 
ine 



An improved knowledge of the emission char- 
acteristics as a function of pulsar rotation phase 
will lead to a better understanding of the emis- 
sion processes and pulsar geometries responsible 
for modulated high-energy 7-rays. In particular, 
the phase-resolved behavior may be used to dis- 
criminate between the two most popular classes 
of models for pulsed high-energy 7-radiation: the 



polar cap models (e.g. Daugherty fc Harding 
1994] |l~996l; ^turner, Dermer, fc Michel 1995[ ; 



Sturner 1995 ), which propose that 7-rays re- 
sult from energetic charged particles accelerated 
along the magnetic field lines just above the po- 
lar cap surfaces, and the outer gap models (e.g. 
Cheng, Ho & Rudcrman 1986a , |b|; Romani &| 



gi daLa hum buLh COS B (pe ar eL al. li^Tj. Yadigaroglu 1995| ; |Yadigaroglu & Romani 1995| ) 



Premer, Hermsen, fc Clear 1^)88) and BGRB1 



([Nolan et al. 1993| ; |Kanbach et al. 1994| ; [Mayer- 
|Hasselwander et al. 1994 ). 



The earlier analyses for Crab and Vela com- 
pared the angular distribution of the data with 
the instrument point-spread distribution to de- 
fine the source flux; no use was made of a de- 
tailed background model. Thus the non-source 
background level was not well determined. In 
the Geminga analysis by Mayer-Hasselwander et 
al. (1994), the non-source background in the 
pulse profile was determined by spatial (cross 
correlation) analysis of the counts in the off-pulse 
phase interval, simultaneously accounting for the 
model-predicted diffuse emission. In the present 
analysis, the increase in photon statistics avail- 
able through the first 3.5 years of EGRET ob- 
servations allows 20 phase intervals to be ana- 
lyzed individually for all three pulsars by a spa- 
tial analysis procedure based on the maximum- 
likelihood concept ( Mattox et al. 1996| ), where 



the non-source background has been accounted 
for on the basis of an improved model of the dif- 
fuse emission ( [Hunter et al. 1997] ). Some differ- 
ences from earlier results may be caused by dif- 
ferences in the analysis techniques. Within this 
analysis each pulsar is analyzed by the same pro- 
cedure, making their results directly comparable. 



in which it is postulated that 7-rays are produced 
by charged particles accelerated in vacuum gaps 
formed in the outer regions of the pulsar magne- 
tosphere. Polar cap 7-radiation should be nar- 
rowly beamed in the same direction as the ra- 
dio emission, with a possible high-energy cutoff 
in the photon spectrum due to attenuation by 
the strong magnetic field near the polar cap sur- 
face. The outer gap 7-ray beam is expected to be 
broad and complex, consisting of emission gener- 
ated along the last closed field lines. 

This paper examines the phase-resolved emis- 
sion characteristics of the three brightest high- 
energy 7-ray pulsars Crab, Vela, and Geminga. 
A limited phase-resolved analysis of the high- 
energy 7-ray emission from PSR B1706— 44 is 
presented by Thompson et al. (1996). The re- 
maining two pulsars detected by EGRET, PSR 
B1055-52 and PSR B1951+32, do not yet have 
sufficient photon counting statistics to perform 
a meaningful phase-resolved study of their emis- 
sion ( [Fierro 1995| ). 

2. Observations 

EGRET is sensitive to 7-rays in the energy 
range from ~ 30 MeV to ~ 30 GeV. The EGRET 
instrument and its calibration have been de- 
scribed extensively in Hughes et al. (1980), Kan- 
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bach et al. (1988, 1989), Nolan et al. (1992), and 
Thompson et al. (1993). EGRET records each 7- 
ray as an electron-positron pair production event. 
This event is processed automatically to deter- 
mine the optimal estimate of arrival direction 
and energy of the photon QBertsch et al. 1989[ ). 
The arrival time of each photon is recorded in 
Universal Coordinated Time (UTC) with an ab- 



umn of Tables |T| and || is the amount of time 
after h UTC on the epoch day that a main ra- 
dio pulse would arrive at the center of the Earth 
if there were no dispersion delays. This quantity 
is used to determine (/>(Tq) in equation (|l|). Since 
Geminga is not a radio pulsar, its timing solution 
is derived using a "Z 2 timing method" based on 



the first four years of EGRET observations (Mat 



sol ite timing accuracy of better than 100 /is. Be- tox 1994 , 1995| ) . The frequency second derivative 



cause of the very low flux level of high-energy 7- 
rays, observing periods typically last about two 
weeks. 

This paper considers the combined observa- 
tions made during the first three cycles of the 
CGRO mission, which lasted from 1991 April to 
1994 October. As shown by Ramanamurthy et 
al. (1995b), Mayer-Hasselwander et al. (1994), 
and McLaughlin et al. (1996), these pulsars show 
no evidence of strong variability. The present 
data include the results of the in-flight calibra- 



tion analysis of EGRET ( Esposito et al. 1997 ). 
This analysis added a small energy-dependent 
correction after the data set used by Ramana- 
murthy et al. (1995b) was complete. 

The photon arrival times were transformed to 
solar system barycentric time T, and the corre- 
sponding rotational phases <j) for each pulsar were 
determined by taking the fractional part of the 
Taylor expansion 

</>(T) = 4>(T ) +fo(T- T ) + \ h (T - T ) 2 

+ |/ 2 (T-T ) 3 , (1) 

where /o , /1 , and fi are the pulsar spin frequency 
and first two time derivatives measured at the 
reference epoch To. The ephemerides used for 
the temporal analysis of the Crab, Geminga, and 
Vela pulsars are listed in Tables |], |||, and ||, re- 
spectively. The epoch is given in units of Modi- 
fied Julian Days (MJD = JD - 2400000.5). The 
Crab and Vela timing solutions are made avail- 
able via the regularly updated Princeton Pulsar 
Timing Database and its associated corrections 
for dispersion measure drifts ([A rzoumanian et al.l 
199^). The value of AT ffi listed in the last col 



/2 of Geminga is too small to be measured and 
is set equal to zero. 

3. Phase-Resolved Spatial Analysis 

Analysis of 7-ray point sources in the EGRET 
dataset is typically performed using a maximum 
likelihood technique ( Mattox et al. 1996| ), 
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which a 7-ray excess above the predicted back- 
ground emission is required to be spatially dis- 
tributed as the instrument point spread function 
before being considered as point source emission. 
In the EGRET energy range, the 7-ray back- 
ground is assumed to consist of an isotropic, ex- 
tragalactic component and Galactic diffuse emis- 
sion due primarily to cosmic-ray particles in- 
teracting with matter and fields in the Galaxy 
QBertsch et al. 1993| ; punter et al. 19971 ). A1 " 
though the background model des not provide 
an adequate global fit to the all-sky data, it is an 
acceptable representation of the sky brightness 
in limited regions if the strengths of the galac- 
tic and extragalactic components are treated as 
adjustable parameters. 

The Crab, Vela, and Geminga pulsars have 
sufficient 7-ray counting statistics that it is pos- 
sible to restrict the spatial analysis to photons 
arriving during a specific phase interval, allow- 
ing the absolute emission levels to be determined 
across the pulse profile. 

The likelihood analysis is subject to small 
systematic uncertainties ( Mattox et al. 1996| ), 
including inaccuracies in the assumed diffuse 
model, assumptions made about other sources in 
the field of view, and limitations of the EGRET 
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Table 1 

Radio Timing Parameters of the Crab Pulsar 





To 


fo 


fi 


h 




AT e 


Valid Dates 


(MJD) 




(HT 10 s- 2 ) (10 


-20 g - 




(ms) 


1QQ1 Mar 16- Tim 5 


48371 


9Q Q4Q951 537Q5Q3 


-3 77657 


n 81 8 




8 3 


1991 May 30-Aug 25 


48449 


29.9467067038240 


-3.77575 


1.06 




23.8 


1992 Jun 28-Oct 18 


48857 


29.9334039601096 


-3.77156 


1.1D 




1.0 


1992 Dec 27-1993 Apr 11 


49035 


29.9276050289288 


-3.76972 


1.05 




25.8 


1993 Apr 6-Jul 5 


49128 


29.9245763662884 


-3.76877 


1.01 




15.5 


1993 Nov 5-1994 Jan 27 


49337 


29.9177729429376 


-3.76650 


1.32 




23.0 


1993 Dec 19-1994 Feb 28 


49375 


29.9165363960966 


-3.76612 


0.921 




4.9 


1994 Jun 23-Nov 15 


49598 


29.9092823758179 


-3.76380 


1.30 




26.3 






Table 2 










Timing Parameters for the Geminga Pulsar 










To 


fo 


h 




h 




Valid Dates 


(MJD) 


(s- 1 ) 




( 






1989 Jul 8-1996 Apr 24 


48750 


4.2176690940300 


-1.95206 x 10" 13 


< 6 


x 10 


-24 






Table 3 










Radio Timing Parameters of the Vela Pulsar 










To 


fo 


h 


h 




AT e 


Valid Dates 


(MJD) 


(s- 1 ) 


(10- 11 s- 2 ) (10 


-21 s - 


3 ) 


(ms) 


1991 Feb 1-May 22 


48343 


11.1988875613749 


-1.55793 


1.18 




66.2 


1991 Aug 12-Sep 29 


48504 


11.1987003711301 


-1.56791 


4.17 




22.3 


1991 Oct 21-1992 Jan 16 


48593 


11.1985799537572 


-1.56584 


2.25 




93.3 


1993 Apr 1-Aug 8 


49142 


11.1978384490139 


-1.56166 


0.876 




84.3 


1993 Jul 1-Aug 27 


49197 


11.1977642483505 


-1.56105 


3.88 




63.5 


1993 Sep 16-Nov 20 


49278 


11.1976550192558 


-1.56063 


1.08 




84.8 


1994 Aug 28-Sep 23 


49605 


11.1972260134043 


-1.56636 


0.626 




69.7 
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calibration data. Thompson et al. (1995), for 
example, recommend a 10% systematic uncer- 
tainty be associated with any absolute number, 
and McLaughlin et al. (1996) find that a 6% 
systematic uncertainty is sufficient to justify the 
assumption that the pulsars are not varying. For 
comparisons between different phases of a given 
pulsar, these systematic uncertainties can be ne- 
glected, because all phases are analyzed under 
the same conditions. 

For each of the three bright 7-ray pulsars, the 
pulse profile is divided into twenty equal-width 
phase bins (A0 = 0.05), and likelihood analysis is 
performed on the twenty resulting spatial photon 
counts maps. There are not enough counts to 
perform a full spectral analysis for each of the 20 
phase bins. Instead, a hardness ratio R is defined 
as 

F(> 300 MeV) 
F(100 - 300 MeV) ' 1 ' 

where F(100 - 300 MeV) is the likelihood flux 
measured between 100 and 300 MeV, and F(> 
300 MeV) is the likelihood flux measured above 
300 MeV, both in units of photons cm~ 2 s . 
A source with a differential flux behaving as a 
perfect power law above 100 MeV with a spectral 
index of —2.0 will have a hardness ratio of R ~ 
0.5. 

3.1. The Crab Pulsar 

Analysis of the Crab region is complicated by 
the fact that a solar flare bright enough to be 
detected by EGRET as a strong source occurred 
on 1991 June 11 ( |Kanbach et al. 1993j ), when the 
Sun was in a direction only ~ 4? 4 from that of the 
Crab pulsar. Rather than introduce possible er- 
rors by trying to account for the solar flare emis- 
sion, all photons arriving during the ~8.5 hour 
flare period are excluded from the Crab pulsar 
analysis. 

The results of the phase-resolved spatial like- 
lihood analysis of the Crab pulse profile are 
shown graphically in Figure [fl. The top panel 
shows the traditional phase histogram above 100 



MeV formed by epoch folding all events from the 
first three cycles of EGRET observations arriving 
within an energy-dependent cone of half-angle 
6 67 = 5?85 x (£ 7 /100 MeV)" - 534 about the pul- 
sar position, where £L is the measured photon 
energy. This cone accepts ~ 67% of the pho- 
tons that EGRET detects from a point source 
( [Thompson et al. 1993| ). The phase interval has 
been divided into the eight components defined 
in Table |], indicated by dashed lines in the top 
panel of Figure |l[ These boundary intervals dif- 
fer from those of Clear et al. (1987) and Nolan 
et al. (1993) because the increase in data has al- 
lowed features to be resolved on a finer scale. 
The horizontal dashed line represents the non- 
source background gamma radiation within this 
cone, determined as follows: The off-pulse phase 
region was analyzed using the maximum likeli- 
hood method to calculate the number of pho- 
tons associated with the Crab. 67% of these 
should be in the light curve. The non-source 
background is then the difference between the to- 
tal number of photons and the number of source 
photons expected in this part of the light curve. 
This method is similar to the one used by Mayer- 
Hasselwander et al. (1994) for Geminga. 

The center panel of Figure [l] shows the instan- 
taneous photon flux above 100 MeV measured for 
each of the twenty phase bins. The total level of 
emission from Crab is the average of the plotted 
values. The error bars represent the la statisti- 
cal uncertainty levels. The values are plotted on 
a logarithmic scale because the photon flux varies 
by almost two orders of magnitude across the 
pulse profile, with the single bin centered about 
the first peak accounting for ~ 35% of the to- 
tal Crab emission. Unlike the phase histogram 
shown in the top panel, the center plot shows 
only the emission above the predicted 7-ray back- 
ground level. Significant point-source emission is 
detected from the Crab region over almost the 
entire pulse profile. The low bin, centered on 
phase —0.15, has statistical significance of only 
2.5a, an excess of 27 photons on the expected 
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Table 4 

Crab Pulsar Phase Interval Definitions 



Component 




Name 


Phase Interval 


Leading Wing 


1 


LW1 


0.88 - 0.94 


Peak 1 




PI 


0.94 - 0.04 


Trailing Wing 


1 


TW1 


0.04 - 0.14 


Bridge 




Bridge 


0.14 - 0.25 


Leading Wing 


2 


LW2 


0.25 - 0.32 


Peak 2 




P2 


0.32 - 0.43 


Trailing Wing 


2 


TW2 


0.43 - 0.52 


Offpulse 




OP 


0.52 - 0.88 



non-source emission of 140. All the other bins 
have significance greater than 3.5<r. 

Previous studies of the high-energy 7-ray emis- 



sion from the Crab complex (Clear et al. 1987; 
Nolan et al. 1993] ; Ramanamurthy et al. 19951) ; 
ie Jager et al. 1996| ) have treated the emission 



in the offpulse (OP) region of the pulse profile as 
an unpulsed component due entirely to the Crab 
nebula. This is in analogy with observational 
results at lower energies, where the total Crab 
luminosity is dominated by the nebular contribu- 
tion. Unlike detectors at lower energies, however, 
current 7-ray telescopes do not have sufficient an- 
gular resolution to separate the pulsar radiation 
from any possible nebular emission. As will be 
shown in the following sections, it is very likely 
that much of the emission in the offpulse phase 
interval comes from the pulsar itself. It is note- 
worthy that the photon flux level measured for 
the phase bin at eft ~ 0.2 is comparable to the flux 
level measured over the offpulse phase interval. 
It had been assumed that the emission through- 
out the bridge region was well above the offpulse 
emission level. An unpulsed component present 
in the Crab pulse profile can be no greater than 
flux measured in the weakest phase bin. Thus, 
an unpulsed component has a maximum flux of 
(3.8 ± 1.5) x 10 -7 photons cm' 2 s -1 , or 16 ± 8% 



of the total Crab emission above 100 MeV. 

The lower panel of Figure [l] plots the hard- 
ness ratios calculated for the same set of phase 
bins. There is no data point plotted for the bin 
centered at phase <j> pa —0.15 because the paucity 
of counts in that bin prevented a detection in 
the 100-300 MeV energy range. The most strik- 
ing aspect of this plot is the absence of a strong 
correlation to the features in the pulse profile. 
The hardness ratio increases smoothly as the pul- 
sar phase moves across the first pulse, leveling 
out in the bridge region before reaching its max- 
imum value in the leading wing of the second 
peak ((f) ~ 0.25). A weak detection of this bin 
in the 100-300 MeV energy range leads to the 
relatively large error bar. As the phase increases 
through the second peak, the hardness ratio de- 
creases back to the soft level of the offpulse emis- 
sion. It is interesting that the extreme hardness 
ratios are in the bridge and offpulse phase inter- 
vals, not near the two peaks which dominate the 
pulsar lightcurve. 

3.2. The Geminga Pulsar 

The phase histogram, phase-resolved instan- 
taneous flux, and hardness ratio distribution of 
Geminga are shown in Figure 0. The vertical 
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Fig. 1. — Phase-resolved likelihood analysis of 
the Crab pulse profile. The top panel shows the 
7-ray phase histogram above 100 MeV formed 
from the data accumulated during the first three 
cycles of EGRET observations. The vertical 
dashed lines indicate the definitions of the pulse 
profile components listed in Table ||. The his- 
togram consists of 100 phase bins. The pul- 
sar phase is defined so that the main radio 
peak occurs at zero. The horizontal dashed line 
shows the level of non-source gamma radiation 
expected, based on spatial analysis of the off- 
pulse region. The center panel shows instanta- 
neous photon flux above the diffuse background 
level as derived from the spatial likelihood analy- 
sis of each phase bin for energies greater than 100 
MeV. The lower panel plots the hardness ratios 
as defined by equation @. 



dashed lines in the top panel indicate the defi- 
nitions of the pulse profile components listed in 
Table ||, and the horizontal dashed line shows 
the non-source background emission, calculated 
in the same way as was done for the Crab. Note 
that the dashed line is close to the observed 
points in the phase region leading the first pulse 



(cf. Fig. 2 of Mayer-Hasselwander et al. 1994). 
The plot of the instantaneous flux level above 
100 MeV shows that Geminga has an even higher 
average level of offpulse emission than does the 
Crab. However, unlike the Crab pulsar, Geminga 
has no known nebular association, and it must 
be assumed that the high level of emission seen 
throughout most of the pulse profile is coming 
directly from the pulsar, as noted previously by 
Mayer-Hasselwander et al. (1994). The flux at 
phase (f) ~ —0.55 is ~ 1/3 the emission mea- 
sured from either of the adjacent phase bins, and, 
even considering the large uncertainties, repre- 
sents a significant drop in intensity. This sta- 
tistical significance of this bin is less than 2.5a, 
so this point should conservatively be treated 
as an upper limit. Any unpulsed component to 
the total Geminga emission can not be stronger 
than the flux level of (5.1 ± 1.7) x 10~ 7 pho- 
tons cm~ 2 s _1 measured in this bin. This corre- 
sponds to 14% ± 7% of the total Geminga emis- 
sion. As opposed to the Crab pulsar, the emis- 
sion level over the entire bridge region of the 
Geminga pulse profile is much higher than the 
offpulse emission level. Except for the one low 
bin, all of the phase regions show excesses at a 
statistical significance greater than 5a. 

The distribution of the Geminga hardness ra- 
tios, plotted in the lower panel of Figure |2|, ap- 
pears to be correlated to the basic features of the 
7-ray pulse profile. There is perhaps a local max- 
imum in the hardness ratio at the location of the 
first peak, with even harder emission in the latter 
part of the bridge region and the second peak. In 
general, the Geminga hardness ratios are much 
higher than those measured for the Crab pulsar. 
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Table 5 

Geminga Pulsar Phase Interval Definitions 



Component 




Name 


Phase Interval 


Leading Wing 


1 


LW1 


0.51 - 0.62 


Peak 1 




PI 


0.62 - 0.71 


Trailing Wing 


1 


TW1 


0.71 - 0.82 


Bridge 




Bridge 


0.82 - 0.96 


Leading Wing 


2 


LW2 


0.96 - 0.09 


Peak 2 




P2 


0.09 - 0.22 


Trailing Wing 


2 


TW2 


0.22 - 0.30 


Offpulse 




OP 


0.30 - 0.51 



3.3. The Vela Pulsar 

Figure || shows the phase histogram, phase- 
resolved flux, and hardness ratios for the Vela 
pulsar, with the pulse profile components defined 
in Table || indicated by vertical dashed lines in 
the top panel. The horizontal dashed line shows 
the non-source background emission, calculated 
in the same way as was done for the Crab. The 
increased statistics and complicated morphology 
of the Vela pulse profile allows for the introduc- 
tion of two additional components not defined by 
Kanbach et al. (1994). For lack of better names, 
these have been denoted as "interpulse 1" (IP1) 
and "interpulse 2" (IP2), and they correspond 
to the slight plateau after the first peak and its 
trailing edge, respectively. These are not to be 
confused with radio pulsar astronomers' use of 
the term "interpulse" to denote the weaker pulse 
in a two-pulse profile. 

The instantaneous flux of Vela shown in the 
center panel of Figure |3] exhibits a greater vari- 
ation as a function of pulsar phase than ei- 
ther Crab or Geminga. The measured instanta- 
neous flux above 100 MeV exceeds 4 x 10~ 5 pho- 
tons cm" 2 s" 1 for both peaks, almost 20 times 
the average intensity from Crab, reflective of the 
fact that Vela is the brightest compact object in 



the 7-ray sky. Although seven of the eight phase 
bins in the offpulse region have a measured flux 
above the diffuse background, none of the indi- 
vidual bins in the offpulse region measure strong 
detections, and one of the bins has only an upper 
limit. This contrasts strongly with the Crab and 
Geminga, which show significant emission over 
most of the pulsar phase. A likelihood analysis 
of the entire offpulse phase region 0.64-0.02 finds 
a point-source signal at a 5.6a confidence level 
and an instantaneous flux of (4.4 ± 0.9) x 10~ 7 
photons cm -2 s _1 . Thus, there is evidence that 
Vela has emission beyond the phase region tradi- 
tionally associated with the pulsar, even though 
the data do not justify the claim that the emis- 
sion extends throughout the full rotation. If this 
level of emission is assumed present throughout 
the pulse profile, it represents only 4.4% ±0.9% of 
the total emission, which does not conflict with 
the upper limits to unpulsed emission established 
by Grenier et al. (1988) or Kanbach et al. (1994). 
Even though the flux levels measured in the off- 
pulse regions of the Crab and Vela pulsars are 
comparable, the Vela flux values have larger un- 
certainties because the EGRET exposure to Vela 
was ~ 2/3 that of the Crab over the first three 
cycles of instrument operation. 

The lower panel of Figure || shows the distribu- 
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Table 6 

Vela Pulsar Phase Interval Definitions 



\J 111 VJ VJ11 C 11 lj 




INTp tti p 

1 1 dlllt, 


T-*h atip Tnt"PT*val 

.1 llCXOC 111UC1 VCXl 


Leading Wing 


1 


LW1 


0.02 - 0.10 


Peak 1 




PI 


0.10 - 0.14 


Trailing Wing 


1 


TW1 


0.14 - 0.19 


Interpulse 1 




IP1 


0.19 - 0.26 


Interpulse 2 




IP2 


0.26 - 0.36 


Bridge 




Bridge 


0.36 - 0.45 


Leading Wing 


2 


LW2 


0.45 - 0.50 


Peak 2 




P2 


0.50 - 0.57 


Trailing Wing 


2 


TW2 


0.57 - 0.64 


Offpulse 




OP 


0.64 - 0.02 



tion of hardness ratios across the Vela pulse pro- 
file. The weak level of offpulse emission combined 
with the limited EGRET exposure to Vela pre- 
vented a useful determination of hardness ratios 
over the phase range 0.85-0.05. The phase evo- 
lution of the Vela hardness ratio is reminiscent 
of the behavior of the Crab hardness ratio (Fig- 
ure [j]). As with the Crab pulsar, the Vela emis- 
sion becomes harder as the phase moves through 
the first peak and softer as the phase increases 
past the second peak into the offpulse compo- 
nent. 

4. Phase-Resolved Spectra 



To further investigate the implied spectral 
variation function of rotation phase implied 
by the hardness ratio distributions of the three 
bright 7-ray pulsars, differential photon spectra 
over the EGRET energy range are derived for the 
various pulsar components defined in Tables ||, [||, 
and fi Each component spectrum is derived by 
performing standard likelihood spatial analysis 
( Mattox et al. 1996| ) in ten independent energy 
intervals using only those photons within the cor- 
responding phase range. As noted by Esposito 



et al. (1997), the EGRET in-flight calibration 
checks showed that the original calibration below 
70 MeV was in error. The corrections necessarily 
introduce larger uncertainties in the 30-50 MeV 
and 50-70 MeV data points than found in the 
other energy bins. 

4.1. The Crab Pulsar 

The photon spectrum derived from likelihood 
spatial analysis of the total Crab emission is 
shown in Figure |] for the energy range 30 MeV- 
10 GeV. Likelihood analysis is not performed 
above 10 GeV because the low photon statis- 
tics make it difficult to perform effective maxi- 



mum likelihood parameter estimation (see Mat 



tox et al. 1996; ). As indicated by the solid line in 
Figure |], the Crab differential photon spectrum 
dN/dE is consistent with a power law of index 
a = — 2.12±0.03 over the energy range 100 MeV- 
10 GeV, similar to the results of Nolan et al. 
(1993) and Ramanamurthy et al. (1995b). How- 
ever, the first three energy bins from 30-100 MeV 
deviate significantly from the trend established 
by the higher energy bins. A spectral power law 
fit including all ten energy bins results in a x 2 
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Fig. 2. — Phase-resolved likelihood analysis of 
the Geminga pulse profile. The top panel shows 
the 7-ray phase histogram above 100 MeV, with 
the vertical dashed lines indicating the defini- 
tions of the pulse profile components listed in Ta- 
ble H The histogram consists of 100 phase bins. 
The pulsar phase is relative to the timing epoch. 
The horizontal dashed line shows the level of non- 
source gamma radiation expected, based on spa- 
tial analysis of the off-pulse region The center 
panel shows instantaneous photon flux above the 
diffuse background level for energies greater than 
100 MeV, while the lower panel plots the hard- 
ness ratios. 



Fig. 3. — Phase-resolved likelihood analysis of 
the Vela pulse profile. The top panel shows the 
7-ray phase histogram above 100 MeV, with the 
vertical dashed lines indicating the definitions of 
the pulse profile components listed in Table ^. 
The histogram consists of 100 phase bins. The 
pulsar phase is defined so that the main radio 
peak occurs at zero. The horizontal dashed line 
shows the level of non-source gamma radiation 
expected, based on spatial analysis of the off- 
pulse region The center panel shows instanta- 
neous photon flux above the diffuse background 
level for energies greater than 100 MeV, while 
the lower panel plots the hardness ratios. 
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value of 47.5 with 8 degrees of freedom (DOF), 
implying that the probability that the 10 mea- 
sured values are drawn from a photon spectrum 
behaving as a simple power law over the entire 
EGRET energy range is 1.2 x 10~ 7 . Even the 
restricted fit above 100 MeV measures a x 2 °f 
10.7 for 5 DOF, corresponding to a probability 
of 5.9%. 



Spectral Index = -2.12 ± 0.03 
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Fig. 4. — Photon spectrum of the total emission 
from Crab derived using spatial analysis. The 
power law fit is to the data points between 100 
MeV and 10 GeV, with the dashed line indicat- 
ing the extrapolated power law down to 30 MeV. 
The spectral index refers to the differential pho- 
ton spectrum dN/dE. 

Figure || shows the instantaneous photon spec- 
tra derived from likelihood analysis of the phase 
intervals defined in Table ||. The time-averaged 
spectra are obtained by multiplying the flux val- 
ues by the phase width of each component. For 
each component, a spectral fit is only performed 
over the energy interval for which the data points 
were reasonably consistent with a power law dis- 
tribution. Not surprisingly, component spectra 
reflect many of the features suggested by the 
hardness plot of Figure |]. The softest compo- 
nents in the Crab lightcurve are the second pulse 
trailing wing (TW2), offpulse (OP), and first 
pulse leading wing (LW1) components. The pul- 
sar emission is hardest in the bridge and lead- 
ing wing of the second pulse (LW2). The only 



component to show some evidence of a high en- 
ergy turnover is the first peak (PI). Neverthe- 
less, the absence of pulsed emission at TeV ener- 



gies (e.g. Vacanti et al. 1991) indicates that the 
pulsed emission must have a cutoff between the 
EGRET and TeV energy ranges. Although based 
on more extensive data, the principal character- 
istics of the phase-resolved spectral components 
are similar to those found by Nolan et al. (1993) 
for the pulsed emission alone. 

Much like the total spectrum of Figure ||, the 
first pulse trailing wing (TW1), bridge, LW2, and 
OP components all appear to have a softer spec- 
trum at energies less than 100 MeV which devi- 
ates from the power law established above 100 
MeV. This implied spectral break is especially 
apparent for the bridge and offpulse regions, the 
two weakest components of the Crab pulse pro- 
file. To investigate this behavior, the spectra are 
fit with a double power law of the form 



dN 
~dE 



photons cm~ 2 s^ 1 MeV -1 . (3) 



Figure |6| shows the double power law fit to the 
total Crab emission, and Figure ^ shows simi- 
lar fits to spectra of the TW1, bridge, LW2, and 
OP components. The various parameters of the 
double power law fits are listed in Table [7|. The 
X 2 values listed in the sixth column indicate the 
measured photon spectrum for each component 
in Table is well-modeled by emission from two 
independent power laws over the energy range 30 
MeV-10 GeV. The second and third columns of 
Table [7| show that the soft component is consis- 
tently fit with a spectral index of ot\ <^ —4 and 
a flux at 100 MeV of Fi(100 MeV) ~ 10~ 8 pho- 
tons cm -2 s _1 , while the fourth and fifth columns 
show that the harder high-energy emission varies 
for each component. This strongly suggests the 
presence of an unpulsed ultra-soft component to 
the Crab total emission which dominates the pul- 
sar radiation from the weaker pulse profile com- 
ponents at energies less than ~ 100 MeV. 
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Table 7 

Crab Double Power Law Spectral Fits 



Component 


Power Law 1 


Power Law 2 


x 2 

(6 DOF) 


Fi(100 MeV) a 




F 2 (100 MeV) a 




Total 


0.7 ±0.2 


-4.89 ± 0.66 


2.3 ±0.3 


-2.05 ±0.04 


11.04 


TW1 


0.6 ±0.4 


-5.27 ±1.00 


1.3 ±0.4 


-1.86 ±0.13 


6.06 


Bridge 


1.2 ±0.3 


-4.33 ±0.47 


0.1 ±0.1 


-1.42 ±0.36 


10.26 


LW2 


1.6 ±0.4 


-3.97 ±0.69 


0.7 ±0.4 


-1.51 ±0.16 


4.02 


OP 


1.0 ±0.2 


-4.13 ±0.46 


0.2 ±0.2 


-1.76 ±0.29 


2.40 



a Modeled instantaneous differential flux at 100 MeV in units of 10 8 pho- 
tons cm -2 s _1 MeV -1 . 



4.2. The Geminga Pulsar 

The total Geminga spectrum plotted in Fig- 
ure [8| is consistent with a power law from 30 
MeV to 2 GeV of spectral index -1.42 ± 0.02. 
Besides being significantly harder and having no 
ultra-soft low-energy component, the Geminga 
total spectrum is much different than the Crab 
total spectrum in that it shows a sharp spec- 
tral turnover at high energies. The spectral 
result appears consistent with that of Mayer- 
Hasselwander et al. (1994), indicating that the 
differences in analysis technique have no strong 
influence on the results. 

The individual component spectra of Geminga 
are shown in Figure ^. With the exception of 
perhaps the TW2 and offpulse phase intervals, 
all of the components show some indication of a 
turnover at high energies, typically around a few 
GeV. As is suggested by the distribution of hard- 
ness ratios, the P2 component has a significantly 
flatter spectrum than any of the other phase in- 
tervals. Indeed, it is the hardest spectrum yet 
measured by EGRET. The softest components 
are the OP and LW1 phase intervals, which are 
also the only two components that are not de- 



tected at the highest energies. The LW1 and PI 
components appear to exhibit the earliest and 
strongest spectral turnovers. Again, these results 
are completely compatible with those of Mayer- 
Hasselwander et al. (1994), although the defini- 
tions of phase intervals are somewhat different. 

4.3. The Vela Pulsar 

The Vela total photon spectrum is presented 
in Figure 10. Over the energy range 30-2000 



MeV, the spectrum is consistent with a power 
low of index —1.62 ± 0.01. Like Geminga, Vela 
exhibits a spectral turnover above ~ 1 GeV. The 
shape of the spectrum is similar to that found by 
Kanbach et al. (1994), although the formal errors 
on the spectral index do not overlap. Some of this 
may be attributable to time-dependent correc- 
tions applied to the long-term EGRET database 



( Esposito et al. 1997 ) after the work of Kanbach 
et al. (1994). It may also reflect the fact that 
systematic uncertainties become important for 
a bright source like Vela, and such systematics 
are not shown in the statistical error bars. Al- 
though it may appear that the Vela spectrum 
is inconsistent with the fitted power law at en- 
ergies below 70 MeV, the uncertainties associ- 
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Fig. 5. — Instantaneous photon spectra derived 
from likelihood spatial analysis of the Crab pulse 
profile components defined in Table |j. The power 
law fit is indicated by a solid line with the differ- 
ential spectral index a listed for each component. 
The upper limits are 2a. 
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Fig. 6. — Double power law fit to the total Crab 
emission photon spectrum. The fit has a form 
given by equation (||). The dotted lines indicate 
the two independent power laws. 
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Fig. 7. — Double power law fits to the likelihood- 
derived photon spectra of the TW1, bridge, LW2, 
and OP components of the Crab pulse profile. 
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Fig. 8. — Photon spectrum of the total emis- 
sion from Geminga derived using spatial analysis. 
The power law fit is to the data points between 
30 MeV and 2 GeV, with the dashed line indicat- 
ing the extrapolated power law out to 10 GeV. 
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Fig. 9. — Instantaneous photon spectra derived 
from likelihood spatial analysis of the Geminga 
pulse profile components defined in Table EL 



ated with the flux values measured for the low- 
est two energy bins are large, and the deviations 
are not statistically significant. Nevertheless, the 
OSSE detection of unpulsed emission from the 
Vela synchrotron nebula in the 0.061 to ~4 MeV 
range flde Jager, Harding, fc Strickman 1996| ), in- 
dicates that there must be a turn-up of the total 
Vela spectrum near or below the EGRET energy 
range. The low-energy EGRET points are consis- 
tent with the cutoff of the synchrotron spectrum 
near 40 MeV predicted by de Jager, Harding, & 
Strickman. 

The instantaneous differential photon spectra 
of the components defined in Table || are pre- 
sented in Figure 11. For the OP component, 



some energy bins had to be combined so that like- 
lihood analysis could produce a reasonable flux 
estimate. Except for the OP interval, all of the 
Vela pulsar components show a turnover at ener- 
gies on the order of a few GeV, with the LW1 and 
PI spectra exhibiting the earliest deviations from 
power law behavior. Only the TW2 component 
fails to produce a signal at the highest energies. 
The IP2 phase interval is the hardest component, 
as was suggested by the hardness ratio distribu- 
tion of Figure ||. Although the phase bins do not 
match those of Kanbach et al. (1994) exactly, 
the similarities of the two sets of phase-resolved 
spectra are evident. 

5. Discussion 

5.1. Observational Results 

A systematic, phase-resolved analysis of the 
total photon flux above 100 MeV reveals that 
Crab, Geminga, and Vela all have significant 
emission in large portions of the "offpulse" phase 
intervals. As opposed to studies of the Crab re- 
gion, multiwavelength observations of Geminga 
and Vela give no reason to believe that their off- 
pulse 7-ray emission might be due to nebular 
emission. If the Geminga and Vela pulsars are 
producing high-energy 7-ray emission through- 
out much of the rotation phase, then it is reason- 
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Fig. 10. — Photon spectrum of the total emission 
from Vela derived using spatial analysis. The 
power law fit is to the data points between 30 
MeV and 2 GeV, with the dashed line indicating 
the extrapolated power law out to 10 GeV. 
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Fig. 11. — Instantaneous photon spectra derived 
from likelihood spatial analysis of the Vela pulse 
profile components defined in Table |(| 



able to assume that at least some of the emission 
in the Crab offpulse region is also generated by 
the pulsar itself. It is worth noting that for all 
three pulse profiles, the lowest level of emission 
occurs just before the first 7-ray peak. While 
this result is not statistically significant, it may 
be an early indication that the offpulse emission 
is really composed of the extended tail(s) of the 
major 7-ray peaks. 

Although it is not apparent from the phase 
histograms shown in the top panels of Figures |l|, 

and ^, phase-resolved analysis shows that the 
bridge emission is different for each of the three 
bright pulsars. Midway between its two peaks, 
the Crab bridge emission drops to an intensity 
level which is consistent with the weak emission 
in the offpulse region. The Vela bridge emis- 
sion also varies substantially, but remains at a 
consistently high flux level and has a minimum 
very near the second peak. The emission from 
Geminga is both strong and relatively steady 
throughout the bridge region. 

Despite the reduced photon statistics incurred 
by phase-resolved analysis, it is possible to ac- 
quire spectral information in most of the twenty 
phase bins by employing a hardness ratio such as 
that defined in equation (||). The Crab and Vela 
spectral variation are qualitatively similar in that 
the hardness ratios increase smoothly through 
the first peak, have maximum values near the 
bridge region, and drop off after the second pulse. 
The Vela hardness ratios are typically greater 
than those of the Crab pulsar. Also, Vela be- 
comes dramatically harder after the first pulse, 
an effect not seen in the Crab hardness ratios. In 
contrast to Crab and Vela, the hardness ratio for 
Geminga is strongly correlated with the features 
in its pulse profile. There are two broad plateaus 
corresponding to the first peak and bridge region, 
with an even larger hardness ratio at the loca- 
tion of the second peak. For all three pulsars, 
the softest ratios are in the offpulse region. 

Spectra derived from likelihood analysis of the 
individual pulse profile components confirm the 
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behavior suggested by the distribution of hard- 
ness ratios. Of the three bright 7-ray pulsars, the 
softest emission occurs in the leading wing of the 
first Crab pulse which has a measured spectral 
index of —2.65 =L 0.15, with an almost identical 
spectral index measured for the trailing wing of 
the second pulse. The hardest pulsar component 
is the second peak from Geminga, which has a 
measured spectral index of —1.27 ± 0.03. The 
latter part of the Vela interpulse component is 



almost as hard with a measured spectral index 
of -1.38 ±0.04. 

For all three pulsars, the component spectra 
share the same basic features as the total spectra. 
That is to say, the high-energy spectral turnovers 
which are apparent in the Vela and Geminga to- 
tal spectra also are also evident in the various 
component spectra. Likewise, the Crab profile 
components show the same soft emission at low 
energies that is apparent in the total Crab spec- 



tra jandonlytheleadin^ 



idence of a high-energy turnover in its photon 
spectrum. Interestingly, for both the Vela and 
Geminga pulsars, the leading peak spectrum sim- 
ilarly turns over at a lower energy than the other 
pulsar components. 

The four weakest Crab pulsar components are 
dominated by ultra-soft emission at energies be- 
low ~100 MeV. Double power law fits to the 
photon spectra from these four components show 
that the soft emission appears at roughly the 
same level in each component, with an aver- 
age spectral index of —4.3 ± 0.3 and a flux of 
(1.1±0.1) x 10" 8 photons cm~ 2 s" 1 MeV" 1 at 100 
MeV. This ultra-soft emission is also apparent in 
a double power law fit to the total Crab emission, 
clearly indicating that it is an unpulsed compo- 
nent to the Crab 7-ray emission. It is not unrea- 
sonable to expect ultra-soft nebular emission out 
to energies on the order of ~ 100 MeV. The emis- 
sion at X-ray and soft 7-ray energies is consistent 
with synchrotron radiation produced by relativis- 



tic electrons streaming from the pulsar ( Kennel & 
Co t-oniti 19841) . Since the Crab pulsar is unlikely 



to produce electrons with energies greater than 
~ 10 15 -10 16 eV, the synchrotron spectrum will 
begin to rollover at energies above a few MeV. 
The soft emission observed from 30-100 MeV is 
then seen as the tail end of the synchrotron emis- 
sion from the nebula. Detailed modeling of the 
Crab emission by de Jager et al. (1996) confirms 
this qualitative picture. 

The detection of unpulsed emission from the 



Crab complex at TeV energies ( Weekes et al 



1989) implies that there is emission from the 
nebula beyond the synchrotron radiation. The 
TeV emission is explainable in terms of inverse- 
Compton scattering ( de Jager et al. 1996| ), which 
has also been proposed as a possible explanation 
of the harder-spectrum offpulse emission above 
1 GeV. However, based on the level of the TeV 
emission, the expected inverse-Compton contri- 
bution at energies on the order of ~1 GeV is be- 
low the EGRET flux measurements (Aharonian 



fc Atoyan 1995 ). Considering the non-zero off- 



pulse emission seen from the Geminga and Vela 
pulsars, it is very likely that the bulk of the hard 
power-law component of the Crab offpulse emis- 
sion is generated by the pulsar itself. 

5.2. Theoretical Implications 

Original polar cap models of pulsar 7-ray 



emission (e.g. [Harding 1981| ) suggested that double- 
peaked pulse profiles were produces by pulsars 
whose rotation axis was nearly orthogonal to its 
dipole axis. However, this could not explain 
peak separations much different than 0.5 in ro- 
tation phase or strong bridge emission between 
peaks. Since at least five of the six high-energy 



7-ray pulsars have more than one peak (Ra- 



manamurthy et al. 1995a| ; |Fierro 1995| ; |Thomp 



[son et al. 1996 ), it is unlikely that these are 
all orthogonal rotators whose beams cross the 
line of sight to the Earth. More recent incar- 



nations of the polar cap models (e.g. Dermer & 



Sturner 1994 ; ^turner Sz Dermer 1994 ; Daugh 



crty & Harding 1994) have made use of the fact 
that a pulsar with closely aligned rotation and 
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magnetic axes pointed in the direction of the 
Earth will produce a hollow cone of emission from 
a single polar cap and result in two observed 
peaks, corresponding to when the line of sight 
from the Earth enters and exits the hollow cone. 
Yet phase-resolved analysis shows that all three 
bright 7-ray pulsars have emission across much 
of the rotation phase, which for the single pole 
outer gap model demands that the rotation and 
dipole axes have an aspect angle no greater than 
the radius of the polar cap cone of emission. For 
the narrow beaming angles associated with the 
standard polar cap model (e.g. Sturrock 1971 ), 
this suggests that there is an very large popula- 
tion of young nearby neutron stars with closely 
aligned rotation and dipole axes. 

Daugherty & Harding (1996) refined the sin- 
gle polar cap model by allowing the accelerat- 
ing regions above the polar caps to extend out 
to heights on the order of a neutron star radius, 
significantly increasing the solid angle into which 
the cone of 7-ray emission is beamed. Not only 



do $s this make it more likely that a pulsar beam 
will cross the line of sight from the Earth, but 
it relaxes the requirement that the pulsar rota- 
tion and dipole axes be nearly aligned. Under 
these assumptions, Daugherty & Harding (1996) 
are able to closely reproduce the Vela pulse pro- 
file observed by EGRET, matching the narrow 
peak duty cycles, observed peak separation, and 
enhanced bridge emission. They also predict 
there will be a finite level of emission outside of 
the two peaks due to residual high-altitude cas- 
cades. A natural consequence of this scenario 
is that the bridge emission will be harder than 
the peak emission since it comes from the inte- 
rior of the polar cap where the curvature radi- 
ation is less likely to cross magnetic field lines 
and be reprocessed into softer cascade photons. 
The emission outside of the peaks is due to high- 
altitude curvature radiation from electrons which 
have lost much of their peak energy, so this emis- 
sion is expected to be the softest. It is also ex- 
pected that the photon spectra from the peaks 



will turn over at lower energies than the spectra 
from the bridge regions since the peak emission 
will be more heavily attenuated by the magnetic 
field. It has already been noted that the lead- 
ing peaks of the Crab, Vela, and Geminga pul- 
sars have spectra which appear to turn over at 
lower energies than the other pulsar components. 
While the measured hardness ratio distributions 
of the Crab and Vela pulsars qualitatively agree 
with the behavior expected for this polar cap sce- 
nario, the phase-resolved spectral behavior of the 
Geminga pulsar presents some problems. By far 
the hardest emission from Geminga comes from 
the second peak. Moreover, the emission between 
the peaks shows a sharper high-energy spectral 
turnover than the second peak. Both of these 
conditions indicate that the emission producing 
the second Geminga peak is not being signifi- 
cantly reprocessed to softer photons by the mag- 
netic field. 

The emission pattern of the magnetic Compton- 
induced pair cascade model ( Dermer fc Sturnerj 



1994 ; Sturner fc Dermer 1994 ) is similar to that 



of the curvature radiation-induced pair cascade 
model of Daugherty & Harding (1996), except 
the high-energy 7-ray emission is produced close 
to the surface of the star, resulting in a much nar- 
rower 7-ray beam. In calculating the 7-ray spec- 
tral index as a function of pulsar phase, Sturner, 
Dermer, &: Michel (1995) find that the hardest 
emission should come from the bridge region and 
the softest emission from the two peaks. This 
agrees somewhat with the observational results 
for the Crab and Vela pulsars — although the ob- 
served hardness ratios are not directly correlated 
with the peaks as suggested by Sturner et al. 
(1995) — but is in conflict with the extremely 
hard photon spectrum measured for the second 
7-ray peak of Geminga. Moreover, the pulse pro- 
files modeled by Sturner et al. (1995) do not re- 
flect the narrow duty cycles of the peaks, strong 
bridge emission, nor the significant offpulse emis- 
sion observed from the three brightest 7-ray pul- 
sars. 
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Chiang & Romani (1992) found that outer gap 
emission from a single pole produces a broad, 
irregularly-shaped beam of emission which is par- 
ticularly dense near the edges, so that two 7-ray 
peaks are observed when the line of sight from the 
Earth crosses these enhanced regions of the 7-ray 
beam, while the inner region of the beam pro- 
vides a significant amount of emission between 
the peaks. A wide range of peak phase separa- 
tions can be accommodated with a proper choice 
of the observer co-latitude. Double-peaked pulse 
profiles of varying phase separation and strong 
bridge emission occur naturally in this model. 
The non-zero offpulse emission can be seen as re- 
sulting from residual pair cascades high in the 
outer magnetosphere Due to the complex ge- 
ometry and interaction of emission processes in 
the outer magnetosphere, no detailed spectra or 
luminosities have been modeled for single pole 
outer gap model. 

A clue to the nature of pulsar 7-ray emission 
may come from the fact that for each of the three 
bright 7-ray pulsars, the individual component 
spectra have comparable high-energy turnovers, 
if any. Magnetic field attenuation would be a 
likely cause of the spectral turnovers if the the 
various pulsar components are being generated 
close to the surface of the neutron star, where 
the local magnetic field is strong. However, if 
the turnovers are due to attenuation by the mag- 
netic field, the Crab pulsar with an inferred sur- 
face magnetic field of 3.7 x 10 12 G should show 
a sharper spectral break than a weaker field pul- 
sar like Geminga, which has an inferred surface 
magnetic field of 1.6 x 10 12 G. It is also pos- 
sible that the spectral turnovers are the result 
of a cutoff in the source distribution of charged 
particles at high energies. Assuming the maxi- 
mum energy attainable by a charged particle is 
directly related to the maximum potential drop 
of the pulsar, which goes as A<I> oc BP~ 2 , where 
B is the surface magnetic field and P is the ro- 
tation period of the pulsar, one would expect the 
Crab pulsar to have the highest charged parti- 



cle distribution turnover energy among the three 
bright 7-ray pulsars, while Geminga should have 
the lowest. This could be why the Crab photon 
spectrum shows no strong spectral turnover be- 
low 10 GeV, while Geminga has a sharp turnover. 
On the other hand, PSR B1951+32, with a A$ 
lower than Vela's, shows no high-energy cutoff 
(Ramanamurthy et al. 1995a). 

A spectral turnover caused by a dropoff in the 
distribution of primary charged particles at high 
energies will be reflected in both the total emis- 
sion photon spectrum and the individual com- 
ponent spectra. However, it is not clear why the 
first peak of the the 7-ray pulsars seems to have a 
slightly stronger cutoff than do the other pulsar 
components. Perhaps the first peak originates 
closer to the surface of the pulsar where the local 
magnetic field is stronger and there is a higher 
degree of magnetic field attenuation, but then it 
might be expected that the first peak should be 
softer than the other components due to the sig- 
nificant reprocessing of high-energy photons to 
lower energies. While Vela and Geminga have 
somewhat softer leading peaks, the Crab pulsar 
does not. 

Even though there is some question as to the 
exact location in the magnetosphere where each 
pulse profile component is generated, the signifi- 
cant variation in hardness ratio as a function of 
rotation phase indicates the line of sight to the 
pulsar is sweeping across field lines with a broad 
range in curvature radii. Most likely, the softest 
components come from the regions with the most 
curved field lines. The two peaks of the Crab and 
Vela pulsars are softer than the bridge regions, 
which suggests that the peaks might be at least 
partially due to the increase in curvature radia- 
tion one would expect along the most curved field 
lines. Yet, the second peak of Geminga, which is 
actually the stronger peak, has the hardest mea- 
sured spectral index of any object or pulsar com- 
ponent observed by EGRET. Either the primary 
charged particle distribution responsible for this 
peak has a very enhanced, flat number spectrum, 
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or, more plausibly, the second Geminga peak 
results from inverse-Compton scattering, which 
can reasonably produce an intense, hard photon 
spectrum. The second peak of Geminga appears 
to be direct evidence that processes besides cur- 
vature radiation and synchrotron radiation are 
playing an important role in the pulsar magne- 
tosphere. 
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